I 


ad-4156  993  4  DIFFERENT  141  PREDICTION  METHOD  FOR  THREE-DIMENSIONAL 

L4NIN4R  AMO  TURBUI..IUI  04910  M  TAYLOR  NAVAL  SHIR 
RESEARCH  AND  DEVELOPMENT  CENTER  BET.. 

UNCLASSIFIED  N  C  DROVES  ET  Al .  JUN  85  DTNSRDC- 65/058  F/0  20/4  NL 


AD-A156  993 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 
UNCLASSIFIED 

2a  SECURITY  CLASSIFICATION  AUTHORITY 
2b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 
4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

DTNSRDC-85/ 058 


REPORT  DOCUMENTATION  PAGE 

™"™~  lib  RESTRICTIVE  MARKINGS 


3  DISTRIBUTION/ AVAILABILITY  OF  REPORT 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION 
IS  UNLIMITED 

S  MONITORING  ORGANIZATION  REPORT  NUMBERS) 


6a  NAME  OF  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 
David  W.  Taylor  Naval  Ship  <lf  applicable) 

R&D  Center  Code  1542 

6c  ADDRESS  (City,  State,  and  HP  Code)  7b  ADDRESS  (C/fy.  State,  and  /IP  Code) 


Bethesda,  Maryland  20084-5000 


8a  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 


8b  OFFICE  SYMBOL 
(If  applicable ) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


Be  ADDRESS  (City  Stile  and /IP  Code) 


10  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM  PROJECT  TASK  WORK  UNiT 

ELEMENT  NO  NO  NO  ACCESSION  Nu 

_  61152N  ZR  0230KU  UN  478002 

11  TITLE  (Include  Security  Classification) 

A  DIFFERENTIAL  PREDICTION  METHOD  FOR  THREE-DIMENSIONAL  LAMINAR  AND  TURBULENT  b'tUNDAR'' 
_ LAYERS  OF  ROTATING  PROPELLER  BLADES 

12  PERSONAL  AUTHOR(S)  ”  ~  ~  ~  . —  — —  ~ . 

Nancy  C.  Groves  and  Ming-Shun  Chans _ 

Ma  TYPE  OF  REPORT  Ub  TIME  COVERED  14  DATE  OF  REPORT  (Year,  Month.  Day)  1'  PAGE  COUNT 

Final _  from _ TO _  1985  June _  25 _ 

’6  supplementary  notation 

Presented  at  the  15th  ONR  Symposium,  Hamburg,  Germany,  3-7  September  1984 

I 

'■  7  _  COSATl  CODES _ 18  SU8JFCT  TfRMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

_  T_Ll2. _ group - su8;GRQup  Three-dimensional,  boundary  layer,  propeller  blades, 

_  -  laminar,  turbulent,  numerical  calculation 

"7  ABSTRACT  (Continye^on^cpye/se  if ^pecesiary  and  identify  by  block  number) 

A  genera^formulation  is  given  for  the  three-dimensional  boundary-layer  flow  on  a 
rotating  propeller  blade.  The  basic  equations  are  presented  in  a  nonorthogonal  co¬ 
ordinate  system  which  rotates  with  the  blade.  Finite  difference  methods  are  used  to 
develop  a  computer  code  for  solving  the  laminar  and  turbulent  boundary- layer  equations. 
The  Reynolds  stress  tensor  is  modeled  by  an  algebraic  eddy-viscosity  formulation.  In 
general,  the  equations  are  solved  numerically  using  the  standard  Keller  box  method. 
However,  regions  of  flow  reversal  across  the  boundary-layer  are  computed  by  the 
characteristic  box  method.  A  companion  geometry  computer  code,  developed  to  model 
propeller  geometry  characteristics,  and  an  existing  inviscid  flow  code  for  computing 
propeller  blade  pressures  are  combined  with  the  boundary- layer  computer  code  to  form^, 

(Continued  on  reverse  side) 


61152N 


COSATl  CODES  | 

GROUP 

SU8-GROUP 

20  Distribution  'availability  of  abstract 

^UNCLASSIFIED*  UNLIMITED  □  SAME  AS  RPT 
,','a  NAME  OF  RESPONSIBLE  INDIVIDUAL 


^^Nanc^C^Iroves^ 
OO  FORM  1473,84  mar 


21  ABSTRACT  SECURITY  CLASSIFICATION 

□  otic  USERS _ UNCLASSIFIED _ _ 

22b  TELEPHONE  (Include  Area  Code)  [  22v  OFFICE  SvMBOl 

(202)  227-1605  Code  1542 


83  APR  edition  may  be  used  until  e«hausted 
All  other  edition*  are  obsolete 

b~ 


_SECURlTY  CLASSIFICATION  Or  t'Ai,F 
UNCLASSIFIED 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


(Blopk  19  Continued) 

— ^an  efficient  computation  scheme.  For  a  given  potential-flow  solution,  a. typical 
boundary-layer  solution  of  690  grid  points  requires  64  seconds  CPU  time  on  a 
CYBER  176  computer.  , 

Computed  results  are  presented  for  several,  blade  geometries.  The  rotating 
segment  solution  compares  well  with  analytical  and  experimental  data.  Predictions 
for  a  model  propeller  also  compare  favorably  with  experimental  data  and  illustrate 
that  two-dimensional  theory  may  provide  adequate  predictions  for  applications 
where  crossflow  effects  are  not  important.  Geometry  effects  of  warp  and  skew  are 
shown  to  be  small  for  the  boundary-layer  predictions  on  three  model  propellers. 


ilor  For 


t  Cflt  i  ~,r) 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


TABLE  OF  CONTENTS 


Page 

LIST  OF  FIGURES .  iii 

NOTATION  .  v 

ABSTRACT  .  I 

1.  INTRODUCTION  .  1 

2.  COORDINATE  SYSTEMS  .  2 

3.  BOUNDARY  LAYER  EQUATIONS  .  3 

4.  INITIAL  CONDITIONS  .  5 

5.  TURBULENCE  MODEL  .  6 

6.  NUMERICAL  ANALYSIS  .  6 

7.  COMPUTATIONAL  RESULTS  AND  DISCUSSION  .  8 

7.1  ROTATING  SEGMENT  .  8 

7.2  PROPELLER  4119 .  11 

7.3  SKEWED  AND  WARPED  PROPELLERS  .  13 

8.  CONCLUSIONS .  14 

ACKNOWLEDGEMENTS  .  15 

REFERENCES .  15 

LIST  OF  FIGURES 

1  -  Schematic  of  Blade  Coordinate  Systems  .  1 

2  -  Computational  Procedure  . 

3  -  Geometry  and  Coordinates  of  Rotating 

Helical  Segment  .  8 

4  -  Laminar  Flow  Comparisons  for  Rotating  Helical 

Segment  at  r/R  =0.95 .  9 

5  -  Variation  of  Momentum  Thickness  6^  on 

Rotating  Helical  Blade  .  .  .  9 


iii 


Page 


6  -  Variation  of  Limiting  Streamline  Angle  $  on  Rotating 

Helical  Segment  .  10 

7  -  Angles  of  Surface  Streamlines  on  a  Rotating  Disk .  10 

8  -  Development  of  the  Streamwise  Skin-Friction 

Component .  10 

9  -  Velocity  Profile  Comparison  on  Model 

Propeller  4119 .  11 

10  -  Variation  of  Shape  Factor  on  Model  4119 .  12 

11  -  Variation  of  Displacement  Thickness 

on  Model  4119 .  12 

12  -  Variation  of  Streamwise  Skin  Friction  Coefficient 

on  Model  4119 .  12 

13  -  Variation  of  Streamwise  Momentum  Thickness 

on  Model  4119 .  13 

14  -  Variation  of  Limiting  Streamline  Angle 

on  Model  4119 .  13 

15  -  Geometry  of  Three  Model  Propellers  .  13 

16  -  Variation  of  Streamwise  Skin  Friction  Coefficient 

on  Three  Model  Propellers  .  14 


iv 


/ 


i 


NOTATION 


Cf0  *  Cfr 


Turbulence  modeling  damping  factor 

Streamwise  and  radial  skin  friction  coefficients,  respectively 


f\  g' 


hl’  h2’  h3 


Kl»  K2 


K12’  K21 


through 
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u,  v,  w 
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ref 

x,  y,  z 
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Nondimensional  velocities,  f '  =  —  and  g'  =  - 

u  u 

e  ref 

Shape  factor 

Metric  coefficients  for  the  £,  ri,  and  £  coordinates, 
respectively 

Geodesic  curvatures  of  the  curves  n  =  constant  and 
£  =  constant,  see  Eq.  4 

Geometric  parameters,  see  Eq.  5 
Mixing  length  approximation 

Solution  coefficient  terms,  see  Eqs.  16  and  17 

Static  pressure  in  the  fluid 

Reynolds  number  at  0.7-radius,  see  Eq.  27 

Velocity  components  in  the  £,  £,  and  q  directions,  respectively 

Reynolds  stresses 
Reference  velocity 

Cartesian  coordinate  system  where  x  is  the  axis  of  rotation 
(positive  displacement  measured  aft),  z  is  upward  positive, 
and  y  forms  a  right-handed  orthogonal  system,  see  Figure  1 

Fraction  of  chord,  measured  from  leading  edge 
Nondimensional  radius,  fraction  of  tip  radius 

Angle  between  the  £  and  n  surface  coordinates 
Limiting  streamline  angle 
Boundary  layer  thickness 
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Displacement  thickness 

Eddy  viscosity  factor  in  turbulence  model,  see  Eq.  21 

Directional  cosines  between  the  systems  (x,  y,  z) 
and  (£,  n»  O 

Streamwise  momentum  thickness 
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K,  n,  ? 
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Fluid  kinematic  viscosity 

Blade  surface  coordinate  system,  £  is  the  nondimensional  fraction 
of  chord  measured  from  the  leading  edge,  n  is  the  nondimensional 
radius  as  a  fraction  of  tip  radius,  and  t,  is  the  outward  normal 
to  the  blade  surface,  see  Figure  1 

Fluid  density 

Transformed  boundary  layer  coordinates,  see  Eq.  10 
Shear  stress  in  the  £  and  n  directions,  respectively 


^ ,  <{)  Two-component  vector  potential,  see  Eq.  12 

H  Constant  angular  velocity  of  the  propeller 
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A  DIFFERENTIAL  PREDICTION  METHOD  FOR  THREE-DIMENSIONAL 
LAMINAR  AND  TURBULENT  BOUNDARY  LAYERS  OF 
ROT ATI NC  PROPELLER  BLADES 


Nancy  C.  Groves  and  Ming  S.  Chang 

David  Taylor  Naval  Ship  Research  and  Development  Center 
Bethesda,  Maryland  20084-5000 


ABSTRACT 

A  general  formulation  is  given  for 
the  three-dimensional  boundary-layer  flow  on  a 
rotating  propeller  blade.  The  basic  equations 
are  presented  in  a  nonorthogonal  coordinate 
system  which  rotates  with  the  blade.  Finite 
difference  methods  are  used  to  develop  a  com¬ 
puter  code  for  solving  the  laminar  and  tur¬ 
bulent  boundary-1 ayer  equations.  The  Reynolds 
stress  tensor  is  modeled  bv  an  algebraic  eddy- 
viscosity  formulation.  In  general,  the 
equations  are  solved  numerically  using  the 
standard  Keller  box  method.  However,  regions 
of  flow  reversal  across  the  boundary-1 ayer  are 
computed  by  the  characteristic  box  method.  A 
companion  geometry  computer  code,  developed  to 
model  propeller  geometry  charac teri s t i cs ,  and 
an  existinq  inviscid  flow  code  for  computing 
propeller  blade  pressures  are  combined  with  the 
boundarv-1 ayer  computer  code  to  form  an  effi¬ 
cient  computation  scheme.  For  a  given 
potential-flow  solution,  a  typical 
boundary- 1  ayer  solution  of  6%  grid  points 
requires  64  seconds  CPU  time  on  a  CYBER  17b 

computer. 

Computed  results  are  presented  for  several 
blede  geometries.  The  rotating  segment  solu¬ 
tion  compares  well  with  analytical  and  experi¬ 
mental  data.  Predictions  for  a  model  propeller 
also  compare  favorably  with  experimental 
data  and  illustrate  that  two-dimensional  theory 
may  provide  adequate  predictions  for  applica¬ 
tions  where  crossflow  effects  are  not  impor¬ 
tant.  geometry  effects  of  warp  and  skew  are 
shown  to  he  small  for  the  boundary- 1 ayer  pre- 
dHt'ons  on  three  model  propellers. 

1.  INTRODUCTION 

The  ability  to  predict  local  flow  proper¬ 
ties  on  propeller  blades  will  aid  in  the  basic 
understandi ng  of  propul sor  performance  and 
cavitation.  Knowledge  of  propeller  viscous 
phenomena  at  model  scale,  such  as  laminar  and 
turbulent  separation,  prior  to  testing,  will 
lead  to  selection  of  test  conditions 
appropriate  to  simulate  prototype  performance, 
uetermi nation  of  the  flow-field  velocity 
distributions  near  the  tip  will  improve  the 


basic  understanding  of  propeller  tip  vortex 
cavitation  and  provide  Information  for  defining 
the  position  of  the  shed  vortex  sheet  in  the 
inviscid  propeller  model. 

This  paper  describes  a  general  method  for 
calculating  the  three-dimensional  boundary 
layer  around  propeller  blades.  The  solution  is 
divided  into  three  calculations  (1)  potential 
flow,  (2)  geometric  parameters,  and  ( 3 j  boun¬ 
dary  layer  flow.  The  potential  flow  computer 
code  adopted  is  the  Brockett  (1981)  code  for 
use  at  design  conditions.  The  geometry  and 
boundary  layer  codes  are  modifications  of  the 
codes  developed  by  Cebeci,  et  al .  (1978)  for 
calculating  three-dimensional  laminar  and  tur¬ 
bulent  boundary  layers  on  ship  hulls.  The 
significant  modifications  to  the  Cebeci,  et  al . 
(1978)  formulation,  for  propeller  blade 
boundary-layer  applications,  include  the  addi¬ 
tion  of  the  Coriolis  and  centrifugal  forces  due 
to  rotation,  the  use  of  a  coordinate  system 
appropriately  describing  propeller  surfaces, 
and  the  specification  of  appropriate  initial 
condi tions. 

The  present  boundary-1 ayer  formulation 
predicts  both  laminar  and  turbulent  flow  using 
the  differential  solution  method.  The  earlier 
calculation  methods  of  both  Droves  (1981),  for 
propeller  blades,  and  Arakawa,  et  al .  (1983), 
for  axial  flow  pump  blades,  compute  only  tur¬ 
bulent  flow  using  a  momentum  integral  solution 
method.  There  is  'much  discussion  regarding  the 
merits  of  each  of  the  solution  methods.  The 
momentum  integral  method  requires  considerably 
less  computation  time.  The  differential 
method  is  generally  considered  to  give  a  more 
accurate  and  complete  prediction  of  the  flow 
characteristics.  However,  the  major  advantage 
of  the  present  solution  procedure  is  its  capa 
bility  of  predicting  laminar  flow.  This  region 
is  particularly  important  for  model  propeller 
applications  where  a  large  region  of  the  flow 
over  the  blades  is  laminar.  Additionally, 
laminar  flow  predictions  are  necessary  for 
instabi 1 i ty . 

A  similar  differential  calculation  scheme 
has  been  developed  by  Itoh,  et  al .  (1984)  for 
predicting  the  three-dimensional  laminar  and 
turbulent  boundary-1  ayer  flow  on  the  rotating 
blades  of  axial  flow  pumps.  The  present  calcu- 
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latlon  procedure  varies  from  their  method  in 
two  major  respects.  First,  Itoh,  et  al .  (1984) 
use  an  orthogonal  coordinate  system  to  repre¬ 
sent  an  axial  flow  pump  blade.  While  this 
system  simplifies  the  governing  flow  equations. 
It  encounters  difficulties  in  grid  generation 
for  blades  which  deviate  from  a  fan  shape.  A 
more  general  surface  coordinate  system  defined 
In  terms  of  x*;,  fraction  of  chord  measured  from 
t  »  leading  edge,  and  xp,  fraction  of  tip 
radius.  Is  chosen  for  the  present  study.  The 
use  of  this  coordinate  system  alleviates  any 
problems  associated  with  the  calculation  grid. 
The  second  variation  from  the  work  of  ltoh,  et 
al .  (1984)  Is  the  specification  of  the  initial 
conditions  at  the  leading  edge  and  the  hub. 

Itoh,  et  al .  (1984)  applied  the  usual  leading 
edge  and  hub  conditions  for  three-dimensional 
boundaries  to  compute  the  growth  for  a  rotating 
blade.  Their  computation  diverged.  In  order 
to  obtain  meaningful  results,  they  applied  an 
averaging  process  to  both  the  leading  edge  and 
the  hub  and  developed  an  Iterative  procedure  to 
obtain  these  average  values  within  a  certain 
accuracy.  In  the  present  study,  the  leading 
edge  and  hub  conditions  are  handled  dif¬ 
ferently.  The  methods  proposed  in  this  paper 
dd  not  require  an  Iterative  procedure  and  the 
solution  converges  quickly  for  all  grid  points. 
The  details  of  the  procedures  are  discussed  in 
Section  5. 

The  boundary-layer  equations  are  solved 
numerically  using  the  Keller  (1970)  two-point 
finite-difference  method  and  the  characteristic 
box  procedure  (see  Bradshaw,  et  al .  (1981))  for 
computing  regions  of  reverse  crossflow.  A  tur¬ 
bulence  model  1$  required  for  closure  of  the 
turbulent  boundary-layer  equations.  The  zero- 
order  eddy-viscosl ty  model  with  gradual  tran¬ 
sition,  given  by  Cebeci  and  Smith  (1974),  is 
adopted  here. 

Even  though  the  present  analysis  allows 
computation  of  the  three-dimensional  laminar 
and  turbulent  houndary-layer  on  the  complex  pro¬ 
peller  blade  geometry,  two  desirable  flow  pre¬ 
dictions  are  beyond  the  scope  of  this  study. 
These  are  the  complex  hub  interaction  and  the 
flow  transition  calculation.  No  hub  interac¬ 
tion  is  considered  and  the  location  of  tran¬ 
sition  is  an  input  parameter  to  the  computer 
code  determined  by  the  user.  It  should  be 
noted  that  the  restriction  of  computing  at 
design  conditions  only  is  a  limitation  of  the 
Brockett  (1981)  formulation  and  not  of  the 
boundary-layer  formulation.  The  boundary  layer 
may  be  computed  with  any  potential  flow  and 
issoeiated  blade  offsets.  Finally,  there  is  no 
distinction  in  the  boundary-layer  solution  pro¬ 
cedure  between  the  pressure  and  suction  blade 
surfaces.  The  differences  between  the  surfaces 
are  accounted  for  in  the  potential  flow  and 
geometry  calculations. 

Results  of  the  propeller  blade  boundary- 
layer  calculation  are  presented  for  several 
geometries.  The  first  blade  studied  is  the 
large  chordlength  segment  investigated  experi¬ 
mentally  and  analytically  by  lakshminarayana, 
et  al .  (1972).  The  laminar  predictions  for  the 
blade  are  in  agreement  with  the  computed 
results  of  both  Ranks  and  f,add  (1962)  for  a 


rotating  sector  and  Morris  (1981)  for  a 
rotating  helical  blade.  Overall,  the  turbulent 
predictions  agree  well  with  both  the  experimen¬ 
tal  data  of  Lakshminarayana,  et  al .  (1972)  and 
the  analytical  values  of  Cham  and  Head  (1969) 
for  the  limiting  case  of  a  rotating  disk. 
Discrepancies  occur  in  the  limiting  streamline 
angle  6  which  is  consistently  over- predicted  by 
the  present  method.  The  second  blade  evaluated 
is  a  model  propeller  designed  at  DTNSRDC  by 
Denny  (1968)  and  tested  in  the  DTNSRDC  12- in. 
water  tunnel  by  Jessup,  et  al .  (1984).  Both 
two-  and  three-dimensional  calculations  were 
made  for  this  blade.  The  results  show  that 
two-dimensional  theory  can  adequately  predict 
the  measured  data  as  well  as  three-dimensional 
theory  with  the  exception  of  the  crossflow. 

The  final  geometries  investigated  are  an 
unwarped,  a  warped,  and  a  skewed  model  propeller, 
examined  earlier  by  Groves  (1981),  with  no 
experimental  boundary-layer  data  available. 

The  current  predictions  show  little  variation 
in  the  predicted  local  skin  friction  coef¬ 
ficient  for  the  three  blades. 

2.  COORDINATE  SYSTEMS 

Figure  1  depicts  a  typical  propeller  blade 
and  hub  configuration.  In  practice,  there  are 
N  symmetrical 1y- spaced  identical  blades 
attached  to  the  hub,  but,  for  clarity,  only  one 
is  shown  in  Figure  1.  The  propeller  rotates 
with  the  constant  angular  velocity  u  about  the 
x-axis.  The  Cartesian  x,  y,  z,  coordinate 
system  is  fixed  in  space  and  does  not  rotate 
with  the  blade.  In  this  system,  x  is  the  axis 
of  rotation,  with  positive  displacement 
measured  aft;  z  is  taken  as  upward  positive; 
and  y  forms  a  rj ght- handed  orthogonal  system  as 
shown  in  Figure  1.  In  the  present  calcu1 ition 
scheme,  the  blade  geometry  and  potential  velo¬ 
cities  are  initially  specified  in  the  Cartesian 
system  by  the  Brockett  (1981)  lifting  surface 
computer  code. 


Fig  1  Schematic  of  Blade  Coordtnatr  Systems 
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Although  the  Cartesian  coordinate  system 
could  be  used  for  the  boundary- layer  solution 
procedure,  a  system  fixed  to  the  blade  surface 
and  rotating  with  the  blade  Is  chosen  for  its 
convenient  representation  of  the  complex 
geometry.  The  surface  coordinates  4,  n,  and  4 
are  defined  with  4  equal  to  the  nondlmenslona! 
fraction  of  chord  measured  from  the  leading 
edge;  n  equal  to  the  nondlmenslonal  radius  as  a 
fraction  of  the  tip  radius;  and  4  equal  to  the 
outward  normal  to  the  blade  surface.  These 
coordinates  are  also  illustrated  In  Figure  1, 
where  a  Is  shown  as  the  angle  between  the  4  and 
n  surface  coordinates.  The  velocity  components 
u,  v,  and  w  are  defined  in  the  4,  4,  and  n 
directions,  respectively. 

The  4  and  n  coordinates,  the  Cartesian 
offsets  on  the  blade  reference  surface  ( 1 . e . . 
no  blade  thickness),  and  the  velocity  com¬ 
ponents  In  a  helical  reference  frame  are  stan¬ 
dard  quantities  produced  by  the  Brockett 
lifting  surface  (1981)  computer  code.  In  pre¬ 
paration  for  Its  use  In  the  present  computation 
system,  the  Brockett  code  (1981)  has  been 
modified  to  calculate  the  Cartesian  offsets  on 
the  actual  blade  surfaces  and  the  velocities  In 
the  Cartesian  reference  frame.  The  surface 
offsets  are  needed  to  define  certain  derivative 
quantities  i . e - ,  the  metric  coefficients, 
geodesic  curvatures,  etc.,  used  in  the 
boundary-layer  solution.  Metric  coefficients, 
denoted  by  hi ,  h? ,  and  h3  for  the  4 ,  n ,  and  c 
coordinates,  respectively,  correlate  surface 
distances  between  the  Cartesian  and  surface 
coordinate  systems.  As  is  typical  practice  In 
boundary- layer  formulations,  the  boundary  layer 
thickness  Is  assumed  to  be  small  and  the  4 
coordinate  is  assumed  to  measure  distance  along 
the  surface  normal.  Therefore,  the  metric  113 
Is  set  to  unity  with  no  loss  of  generality. 

The  remaining  metric  coefficients  are  defined 
as 
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In  these  equations,  u,  v,  and  w  represent 
the  velocity  components  In  the  4,  4,  and  0 
directions,  respectively,  a  is  the  angle  between 
the  surface  coordinates  4  and  n,  and  n  Is  the 
constant  angular  velocity.  The  metric  coef¬ 
ficients  hi  and  hj  were  defined  in  the  previous 
section.  The  quantities  Ki  and  K2  are  the 
geodesic  curvatures  of  the  curves  ^constant 
and  4-constant,  respectively,  and  are  defined 
by 
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The  parameters  and  Kp  1  arp  defined  in  terms 
of  the  geometry  as 
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3^  BOUNDARY  LAYER  EQUATIONS 

The  steady,  incompress’hle  boundary 
layer  equations  for  rotational  flow  in  a  non- 
orthogonal  coordinate  system  are  given  by 
Yamataki  (1981).  The  equations  are  identical 
to  those  presented  by  Cebeci,  et  al .  (1478), 
with  tne  addition  of  two  terms  in  the  'nomentur 
equations  representing  the  centrifugal  and 
Coriolis  forces.  The  governing  equations, 
using  the  surface  coordinates  4,  • ,  and  ;,  are: 
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The  remaining  parameters  in  thp  boundary  layer 
equations  are: 

j  :  static  pressure  in  th-  fluid 
,  -  fluid  density 

-  shear  stress  in  the  .-direction. 
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t2  »  shear  stress  In  the  n-dlrectlon,  (6)  where  *>  and  *  are  defined  as 
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v  *  kinematic  viscosity  of  the  fluid 

u'v',  v'w'  *  Reynolds  stresses 
a  *  constant  angular  velocity  of  propeller 
n  ,  n  ,  n  *  directional  cosines  between  the 
123  systems  (x,  y,  z)  and  (£,  n,  cl 

The  boundary  layer  thickness  Is  denoted  by  6 
and  the  boundary  layer  Is  defined  In  the  region 
0<c«6.  The  boundary  conditions  are: 
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The  pressure  gradients  are  determined  from  the 
momentum  equations  at  the  edge  of  the  boundary 
layer.  The  edge  equations  are 
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n-momentum: 
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Transformed  variables  are  defined  as 
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A  two  component  vector  potential  is  introduced 
wi  th 
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and  uref  Is  a  reference  velocity. 

The  potential  velocities  u,  v,  and  w  are 
rewritten  using  the  vector  potentials  a  and  a. 
the  Reynolds  stresses  are  rewritten  In  terms  of 
the  eddy  viscosity  factor  cB,  and  the  boundary- 
layer  edge  conditions  are  applied  to  eliminate 
the  pressure  terms.  The  transformed  e-  and  n- 
equatlons  are 
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In  these  equations,  primes  denote  differen¬ 
tiation  with  respect  to  a  and 
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The  coefficient  terms  inj  through  mi;>,  identical 
to  those  given  by  Cebeci,  et  al .  (1978),  are 
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SUCTION  SURFACE 


Fig.  1C  Variation  of  Straamwiu  Skin  Friction  Coafficiant 
on  THraa  Modal  Propwller* 


chord  value  xj;  =  0.18.  The  design  advance 
coefficient  Jy  of  all  three  models  Is  0,889  and 
the  flow  conditions  are  v  =  1.191  x  10-°  mZ/s 
(1.28  x  10-5  ft’/sec ) ,  V*16.9  m/s  (55.3 
ft/sec),  and  ft  =  391.5  rad/sec. 

figure  16  shows  the  predicted  skin  fric¬ 
tion  coefficient  on  the  suction  surface  for  the 
three  model  propellers  at  the  fraction  of 
radius  value  xp  -  0.91.  This  location 
corresponds  to  streamline  14  used  for  com¬ 
parisons  in  the  earlier  work  of  Groves  (1981) 
where  the  turbulent  boundary  layer  Is  computed 
using  momentum  integral  methods.  The  present 
results  predict  only  a  slight  increase  In  the 
sk'n  friction  coefficient  for  the  warped  and 
skewed  blades.  This  contradicts  the  results 
obtained  by  the  integral  method  of  Groves 
(1981)  where  the  skin  friction  coefficient  was 
predicted  to  Increase  by  nearly  30  percent  for 
the  zero  warp  blade.  An  Investigation  into  the 
momentum  integral  coding  has  Identified  an 
error  in  the  specification  of  certain  geometry 
parameters  for  these  blades  and  thus  leads  to 
the  different  and  erroneous  earlier  conclusion 
of  Groves. 

8.  CONCLUSIONS 

This  paper  presents  analysis  and  results 
for  computing  three-dimensional  laminar  and 
turbulent  boundary  layers  on  the  surface  of  a 
propeller  blade  using  the  differential  method. 
The  solution  procedure  Is  a  modification  of  the 
Cebeci  ,  et  al .  (1978)  scheme  for  boundary- layer 
calculations  on  three-dimensional  ship  hulls. 

The  major  changes  to  the  scheme  result  from  the 
propeller  blade  rotation,  the  complex  propeller 


geometry,  and  the  specification  of  Initial  con¬ 
ditions  In  order  to  begin  the  calculation  pro¬ 
cedure.  Typical  boundary  layer  computations 
use  a  grid  of  30  chordwlse  and  23  radial  points 
thereby  covering  the  entire  blade  surface.  For 
this  grid  size,  64  seconds  CPU  time  on  a  CYBER 
176  computer  are  required  for  a  complete 
calculation. 

Overall,  the  predicted  boundary-layer 
parameters  are  shown  to  give  reasonable 
agreement  with  experimental  data  for  both 
simple  (rotating  disk  and  rotating  helical 
segment)  bodies  and  a  model  propeller.  The 
largest  discrepancies  between  the  current  pre¬ 
dictions  and  experimental  data  and  previous 
theories  occur  for  the  limiting  streamline 
angle  g  which  Is  overpredicted  In  the  turbulent 
region.  For  the  rotating  segment,  this  para¬ 
meter  Is  shown  to  be  correctly  computed  by  the 
differential  method  In  the  laminar  region.  For 
turbulent  flows,  momentum  Integral  methods, 
which  Include  a  modified  entrainment  function 
to  account  for  rotational  effects,  have  been 
shown  (see  Banks  and  Gadd  (1962)  and  Groves 
(1981))  to  Improve  the  prediction  of  tang.  The 
viscosity  coefficient  In  the  present  differen¬ 
tial  method  has  not  been  modified.  The  value 
of  tang  may  be  reduced  if  the  eddy  viscosity 
were  modeled  more  precisely. 

The  analysis  of  model  Propeller  4119  has 
provided  Insight  into  the  boundary-1 ayer  calcu¬ 
lation  on  propeller  blades.  First,  the  large 
extent  of  laminar  flow  on  the  model  blade 
necessitates  the  capability  of  a  laminar  calcu¬ 
lation  procedure  and  an  Instability  prediction 
method  to  determine  the  transition  location. 
Figures  9  through  13  show  that  two-dimensional 
theory  can  be  used  for  the  prediction  of 
boundary-layer  parameters  such  as  chordwise 
velocity  profiles,  shape  factor,  and  chordwise 
skin  friction.  However,  applications  involving 
crossflow,  such  as  the  Instability  calculation 
or  tip  vortex  Investigations,  require  the  full 
three-dimensional  computation.  Figures  9  and 
14  show  the  dangers  of  using  a  partially 
rotating  three-dimensional  calculation. 

Incorrect  values  of  the  crossflow  velocities 
and  tang  are  predicted  when  rotation  is 
included  in  the  potential-flow  solution  and  not 
included  in  the  boundary-layer  solution. 

Calculations  of  the  three-dimensional 
boundary-layer  characteristics  of  three  model 
propellers  of  varying  geometry  indicate  only 
slight  differences  in  the  skin  friction 
parameters;  see  Figure  16.  This  result, 
contrary  to  the  earlier  momentum  integral 
result  of  Groves  (1981),  led  to  the  discovery 
of  a  geometry  error  in  the  earlier  computer 
code. 

The  results  in  this  paper  are  encouraging. 
Although  computed  with  the  present  preliminary 
version  of  the  code,  comparisons  with  measured 
data  are  quite  good.  Further  code  modifica¬ 
tions  are  anticipated  to  allow  the  calculation 
to  proceed  past  the  location  of  turbulent 
separation  and  to  improve  the  eddy  viscosity 
model.  Additional  comparisons  will  be  made 
with  the  experimental  data  of  Itoh,  et  al  . 

(1984)  for  a  rotating  axial  flow  pump  blade. 
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meter  is  positive  for  flow  with  rotation  and 
negative  for  flow  with  no  rotation,  as  pre¬ 
dicted  for  the  radial  velocity  distribution. 

To  sumnarize,  the  boundary  layer  charac¬ 
teristics  of  a  propeller  blade  can  be  computed 
with  reasonable  accuracy.  Two  dimensional 
theory  gives  excellent  predictions  of  the  flow 
for  applications  where  crossflow  is  not  impor¬ 
tant.  Applications  where  the  crossflow  is 
important,  including  model  propeller  studies 
where  large  regions  of  laminar  flow  exist  and 
in  the  study  of  instability,  require  use  of  the 
full  three-dimensional  equations.  The  use  of 
three-dimensional  theory  without  including 
rotation  in  the  houndary-1 ayer  solution  should 
be  avoided. 


In  this  section,  results  are  presented  for 
the  computed  boundary-layer  characteristics  of 
three  analytically-defined  propeller  blades. 
These  three  0.30  m  (1  ft)  diameter  model  pro¬ 
pellers  were  chosen  to  investigate  the  effects 
of  varying  geometry  on  boundary- layer  flow. 

The  propellers,  denoted  as  Model  4381,  an 
unwarped  blade;  Model  4498,  a  72-degree  warped 
blade;  and  Model  4383,  a  72-degree  skewed 
blade,  are  depicted  in  Figure  15.  The  complete 
geometry  of  these  mode!  propellers  is  given  in 
Groves  (1981)  and  is  not  repeated  here. 

Since  no  experimental  data  exist  for 
these  particular  models,  the  transition  loca¬ 
tion  must  be  estimated.  Unpublished  experimen¬ 
tal  results  at  DTKSRDC  on  a  similar  geometry 
indicate  that  the  flow  over  the  blades  is  fully 
turbulent  at  a  0.7-radius  Reynolds  number  of 
4  x  10®.  This  Reynolds  number  corresponds  to 
tripping  the  boundary  layer  at  the  fraction  of 


la)  Model  4381  with  no  warp 


lb!  Model  4488  with  72  degree  warp 


Id  Model  4381  with  72  degree  skew 


Fig.  15  Geometry  of  Three  Model  Propellers 
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«nd  the  chordwise  velocity  u.  In  addition,  the 
momentum  thickness  9\\  and  skin  friction  coef¬ 
ficient  Cf  are  compared  with  two-dimensional 
predictions  from  a  computer  code  developed  by 
Cebeci  (1978).  The  radial  velocity  w  and  the 
limiting  streamline  angle  a  are  the  two  parame¬ 
ters  shown  which  are  unique  to  three- 
dimensional  flow.  No  measurement  has  been  made 
of  these  quantities  as  yet.  A  second  calcula¬ 
tion  was  made  using  the  present  method.  This 
calculation  is  performed  with  the  rotational 
effects  included  in  the  potential-flow  veloci¬ 
ties  but  not  in  the  boundary- layer  computation. 
These  results,  which  approximate  the  rotating 
flow  above  a  solid  ground,  are  also  presented. 

Figures  9  through  14  show  the  boundary-layer 
comparisons  for  Model  4119.  The  series  of 
dashed  and  dotted  lines  in  Figures  9a  through 
9c  represent  Jessup's  (1984)  measured  chordwise 
velocity  profiles  for  the  three  blades  at  three 
chordwise  locations,  x^  »  0.1,  0.4,  0.8, 
respectively.  The  boundary  layer  is  shown  to 
thicken  considerably  as  the  blade  trailing  edge 
is  approached.  The  open  circles  and  squares, 
denoting  the  fully-rotational  three-dimensional 
calculation  and  the  two-dimensional  calcula¬ 
tion,  respectively,  are  both  shown  to  approxi¬ 
mate  the  measured  chordwise  profiles  equally  as 
well.  However,  the  calculation  without  the 
rotation  in  the  boundary  layer  overpredicts  the 
velocity.  It  is  further  seen  from  Figure  9 
that  disregarding  the  rotational  effects  in  the 
bounlary-layer  computation  leads  to  the  predic¬ 
tion  of  an  Inward  rather  than  an  outward  radial 
flow.  An  inward  flow  contradicts  the  flow 
visual ization  results  of  Jessup  (1984). 
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MEASURED 

8IAOE  1  _ 

BIAOE  2 _ _ 

BLADE  3 _ 

COMPUTEO 
0  3DW  44  RAO  SEC 
:  2  D 

m  3  D  W  0  RAD  SEC 

Fiy.  10  Variation  of  SHapa  Facto*  on  Modal  41  Ik 
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Fig.  12  ■  Variation  of  Straamwiaa  Skin  Friction  Coefficient 
on  Modal  4119 


Considerable  scatter  is  noted  in  Figure  10 
for  the  measured  shape  factor  parameter  H. 
Nevertheless,  the  overall  character  of  the 
variable  is  reasonably  wel 1 -predicted  by  the 
various  theoretical  methods,  including  the 
three-dimensional  boundary  layer  without  rota¬ 
tion.  The  displacement  thickness,  shown  in 
r1gure  11  is  again  well -predicted  by  two- 
dimensional  flow  theory  and  three-dimensional 
flow  theory  with  rotation.  Without  the  rota¬ 
tional  effects,  three-dimensional  flow  is  shown 
to  underpredict  the  displacement  thickness. 

Jessup  (1984)  does  not  present  measure¬ 
ments  for  either  the  streamwise  skin  friction 
coefficient  Cf  or  the  momentum  thickness  tljj. 
Comparisons  of  the  calculated  values  of  Cf  and 
'll  are  given  in  Figures  12  and  13.  Again, 
two-dimensional  theory  agrees  well  with  three- 
dimensional  predictions  with  rotation  and  both 
disagree  somewhat  with  three-dimensional  pre¬ 
dictions  with  no  rotation. 

Finally,  the  tangent  of  the  limiting 
streamline  angle  e  is  compared  in  Figure  14  for 
the  three-dimensional  calculations  with  and 
without  rotation.  Unsurprisingly,  this  para- 
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8,  taken  from  Chen  end  Heed  (1969),  shows  cote- 
perl  sons  of  experimental  data  with  several 
calculation  Methods.  The  circled  asterisks 
denoting  the  present  predictions  are  In  good 
agreeMent  with  both  previous  theories  and 
measured  data. 

In  suamiary,  the  present  differential 
method  solution  procedure  has  been  shown  to 
accurately  predict  both  laminar  and  turbulent 
boundary-layer  characteristics  on  a  simple 
three-dimensional  blade.  The  previous  momentum 
Integral  methods  of  Groves  (1981)  and  Arakawa, 
et  al .  (1983)  could  not  predict  laminar  flow,  a 
region  Important  In  model  propeller  and  Insta¬ 
bility  applications.  The  limiting  streamline 
angle  for  turbulent  flow  Is  overpredicted, 
perhaps  Indicating  that  a  change  In  the  eddy 
viscosity  model  Is  needed  to  account  for  the 
rotation. 

7,2  Propeller  4119 


DTNSRDC  Model  Propeller  4119  (see  Denny 
(1968)  for  the  complete  model  geometry)  was 
evaluated  experimentally  on  the  suction  surface 
by  Jessup,  et  al .  (1984)  in  the  OTNSROC  24-1n 
water  tunnel.  The  three-bladed  unskewed,  pro¬ 
peller  model  has  a  0.30  m  (1  ft)  diameter  with 
a  hub  radius  of  0.03  m  (0.1  ft).  The  design 
advance  coefficient  Jv  of  the  model  Is  0.833 
where  Jv  -  V/(nO)  and  V  equals  the  onset 
speed,  n  Is  the  constant  rotational  speed  In 
revolutions  per  unit  time,  and  0  Is  the  rotor 
diameter.  The  flow  conditions  set  for  Jessup's 
(1984)  experiments  are  the  kinematic  viscosity 
of  68°  fresh  water  y  *  1.00  x  IO-GhK/s  (1.08  x 
10'5  ftz/sec),  the  onset  velocity  V  *  1.83  m/s 
(6.0  ft/sec)  and  the  rotational  speed  n  =  44 
rad/s  (7  rev/sec).  These  conditions  yield  a 
0.7- radius  Reynolds  number  of  Rn  *  7.3  x  10’ 
where 

R„  =  (c)0-7  v 


and  ( c )o . 7  3  0.14  m  (0.46  ft)  Is  the  blade 
chord  at  0.7  radius. 

The  0.7  radius  was  chosen  for  comparison 
with  Jessup's  (1984)  experiments  since  the 
measured  data  are  most  complete  at  this  radius. 
How  visualization  techniques  predict  fully 
turbulent  flow  begins  In  the  region  between 
fraction  of  chord  values  x<-  of  0.5  to  0.6.  The 
present  theoretical  eddy  viscosity  model  Incor¬ 
porates  an  Intermlttancy  region  of  gradual 
transition  from  laminar  to  turbulent  flow.  The 
use  of  an  Intermlttancy  region  requires  an 
early  specification  of  the  transition  location 
to  the  computational  scheme.  It  has  been 
determined  by  trial  and  error  that  an  Input 
transition  location  to  the  computational  proce¬ 
dure  of  X£  *  0.43  yields  fully  turbulent  flow 
at  the  chord  values  of  0.5  to  0.6. 

The  boundary-layer  characteristics  com¬ 
pared  with  experimental  data  are  the  streamwlse 
displacement  thickness 

5*  *  "  1  -  —  dt,  the  shape  factor  h  =  «*/e, 
6  ue 


0  0  2  0  4  0  6  0  8  1  0 

CHORDWISE  VELOCITY  COMPONENT,  u/Ue 
OR  SPANWISE  VELOCITY  COMPONENT.  w/Ue 


OR  SPANWISE  VELOCITY  COMPONENT.  w'Ue 


CHORDWISE  VELOCITY  COMPONENT.  u/U* 
OR  SPANWISE  VELOCITY  COMPONENT.  w/U, 


Fig  9  Velocity  Profile-  im  Mortwl  Pr  ope  Wei  4119 
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the  flow  parameters  calculated  by 
lakshminarayana,  et  al .  (1972). 

figures  5a  and  5b  show  the  momentum 
thickness  as  a  function  of  chordwlse  position 
for  fraction  of  radius,  hr,  values  of  0.72  and 
0.93,  respectively.  The  streamwlse  momentum 
thickness  Is  defined  as 


The  dotted  lines  denoting  the  present  calcula¬ 
tion  and  the  dashed  lines  denoting  the  calcula¬ 
tion  of  Lakshminarayana,  et  al .  are  In  close 
agreement.  The  experimental  data  are  denoted 
by  the  circular  symbols.  At  the  radial  loca¬ 
tion  r/R  =  0.72  shown  In  Figure  5a,  all  calcu¬ 
lation  methods  overpredict  the  measured  values 
of  9ji  for  2  radians  <  0  <  5  radians.  However, 
agreement  between  experiment  and  theory  Is 
quite  good  for  the  early  stages  of  turbulence, 
9x2  radians,  and  for  the  blade  trailing  edge, 
9  •>  5  radians.  The  experimental  data  shown  In 
Figure  5b  at  the  radius  r/R  =  0.93  are  more 
scattered.  Agreement  between  experiment  and 
theory  at  this  radius  is  reasonable  only  for 
9  x  2  radians. 

The  tangent  of  the  limiting  streamline 

angle  p,  j  =  tan-1  f—i, is  shown  in  Figure  6. 

\Cfe/ 
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Fig.  7  .  Angles  of  Surface  Streamlines  on  a  Rotating  Dish. 

X  .  experiment;  A,  Cham  and  Head  calculation  using  the 
entrainment  method;  B,  calculation  with  isotropic  eddy 
viscosity;  C,  Banks  Er  Gadd  (1982);  D,  Goldstein  (1936).  E,  von 
Karman  119211. 


Component. - A.  Cham  and  Head  calculation  using  the 

entrainment  method;— — B,  von  Karman  (19211;  - . 

Goldstein  ( 19351; - ,  calculation,  with  isotropic  eddy 

viscosity;  t  ,  Theodorsen  Er  Regier  119441  Results  taken 
from  Clauser  plots  of  the  present  measurements  x  .516 
rev/m  in;  O  ,  1000  rev/min;  ,  1550  rev/min 


Fig.  9  -  Variation  of  Limiting  Streamline  Angle  (1  on  Rotating 
Helical  Segment 


The  momentum  integral  methods  of  both 
Lakshminarayana,  et  al .  (1972)  and  Droves 
(1981)  show  closer  agreement  to  the  measured 
values  of  tane  than  does  the  present  method. 
Figure  7,  from  the  Cham  and  Head  (1969)  study 
of  the  turbulent  boundary- layer  flow  on  a  disk, 
shows  the  discrepancy  in  various  calculation 
methods  of  the  parameter  tana.  As  shown  by  the 
dotted  line  labeled  R  and  the  circular  symbols 
denoting  the  present  calculation  method,  the 
eddy  viscosity  model  predicts  a  larger  value  of 
the  limiting  streamline  than  do  the  other  tur¬ 
bulence  models.  Referring  to  Figure  6  again, 
the  solid  line  labeled  tana  =  1.8"  shows  the 
exact  analysis  result  of  Ranks  and  Dadd  (13b?) 
for  laminar  flow  on  a  rotating  segment.  The 
present  calculation  agrees  well  with  this  line. 
One  further  factor  to  note  Is  that  the  momentum 
integral  method  of  Droves  (1981)  included  a 
modification  to  reduce  the  entrainment  function 
from  the  flat  plate  value  to  include  the  rota¬ 
tional  effect.  The  viscosity  coefficient  in 
the  differential  method  has  not  been  modified. 
The  current  tang  prediction  may  be  reduced  if  a 
more  precise  eddy  viscosity  model  were  used. 

The  final  comparison  shows  the  skin  fric¬ 
tion  coefficient  for  a  rotating  disk.  Figure 
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(a)  STREAMWISE  SKIN  FRICTION  COEFFICIENT 


(b>  RADIAL  SKIN  FRICTION  COEFFICIENT 
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U.l  VELOCITY  PROFILES 

Ficj  4  L.immar  Flow  Comparisons  for  Rotating  Helical 
Segment  at  r  R  0  95 

Crosswise  and  radial  velocity  profiles  are 
compared  in  Figure  4c.  Morris's  profiles  are 
shown  at  chordwise  locations  of  0,  0.26‘T  and 
0.5?  radians.  The  present  predictions  are 
shown  at  h  :  0.3?  radians  as  the  circle.  The 
square  symbols  denote  the  radial  and  crosswise 
velocity  of  Banks  and  fiadd  (l%?)  at  •»  =  0.1? 
radians.  This  figure  shows  that  the  laminar 
;rg sswise  velocity  profile  is  predicted  quite 
well  and  the  radial  velocity  profile  is  just 
slightly  underpredi c ted  by  the  present  method. 


COMPUTED 


FLAT  PLATE 

LAKSHMINARAYANA  ET  AL 

INTEGRAL 

DIFFERENTIAL 


Fig.  5  Variation  of  Momtntum  Thicknasa  Hn  on  Rotating 
Halical  Blada 


4  second  calculation  was  sade  witn  tran¬ 
sition  set  to  u  -  0.3?  radians.  This  location 
is  forward  of  the  experimental ly  determined 
transition  location  of  1  -  J.M  radians. 
However,  with  the  gradual  transition  model  used 
in  the  current  computer  code,  this  forward 
transition  location  gives  good  agreement  with 
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visualization  techniques,  or  by  assuming  tran¬ 
sition  empirically.  One  empirical  formulation 
which  may  be  used  is  the  flat-plate  critical 
Reynolds  number,  Rx  »  £lr28/v  =  3  x  105.  The 
input  transition  location  is  overridden  if 
laminar  separation,  identified  as  a  negative 
streamwise  skin  friction  coefficient,  occurs. 
Once  transition  occurs,  each  succeeding  node  is 
computed  as  turbulent.  Presently,  the  calcula¬ 
tions  stop  once  turbulent  separation  is  encoun¬ 
tered. 

7.  COMPUTATIONAL  RESULTS  AND  DISCUSSION 

A  variety  of  geometries  were  considered 
during  the  verification  stage  of  the  three- 
dimensional  boundary-layer  computation  scheme. 
These  geometries  include  a  flat  plate,  a 
rotating  disk,  rotation  above  the  ground,  a 
skewed  and  an  unskewed  rotating  segment,  an 
upwarped  mode!  propeller,  and  three  model  pro¬ 
pellers,  one  with  warp,  one  with  skew,  and  one 
with  neither  warp  nor  skew.  Overall  agreement 
with  exact  solutions  and  test  results  is  quite 
good.  The  largest  discrepancy  with  experimental 
data  occurs  for  the  limiting  streamline  angle  o 
in  turbulent  flows. 

Computat i ona 1  results  are  presented  for 
the  following  geometries.  A  single,  nonlifting 
blade  of  large  chordlength  is  examined  ini¬ 
tially.  Predictions  in  the  initial  laminar 
flow  region  are  compared  with  the  results  of 
Banks  and  Gadd  (1962)  on  a  rotating  sector  and 
with  Morris's  (1981)  predictions  on  a  helical 
blade.  The  boundary  layer  on  this  blade  is 
computed  a  second  time  with  a  small  laminar 
flow  region  to  simulate  the  test  conditions  of 
Lak shmi narayana ,  et  al.  (1972).  At  large  chord 
length,  this  blade  approximates  a  rotating  disk 
and  comparisons  are  shown  with  turbulent  flow 
solutions  for  a  disk.  To  investigate  the  lami¬ 
nar  flow  region  more  extensively,  the  flow 
about  a  three-bUded  mode!  propeller  is  exa¬ 
mined.  Transition  tor  this  blade  is  set  at  43 
percent  of  the  chord  length  to  simulate  the 
test  data  of  Jessup,  et  al .  (1984). 

Comparisons  are  made  with  the  suction  surface 
measurements  of  Jessup,  et  al.  at  the  fraction 
o*  radius  xq  =  0.7.  Fina'ty,  the  effects  on 
the  boundary  layer  of  the  geometry  parameters 
warp  and  Skew  are  examined  using  three  model 
propellers.  No  experimental  results  exist  for 
tnese  blades. 

_2_.  I  Pointing  Segment 

The  large  chord  length  rotating  segment 
examined  is  the  blade  generated  for 
investigation  of  turbulent  flow  by 
l.  ak  shmi  narayana ,  et  al.  (1972).  The  blade, 
shown  in  Figure  3,  is  a  single  nonlifting 
rotating  sector  of  92.6  cn  (36.6  in.)  diameter 
with  a  300  degree  included  angle.  The  hub 
radius  is  22.86  cm  (0.76  ft.).  The 
pi tch/di ameter  ratio  of  0.273  for  this  blade  is 
small  enough  to  approximate  the  blade  by  a  flat 
circular  disk  with  a  leading  and  trailing  edge 
300  degrees  apart.  Present  calc  ilat ions  are 
performed  using  the  surface  coordinates  .  (or 
xc,  fraction  of  chord),  ■<  (or  xp  ,  fraction  of 


Fig.  3  -  Geometry  end  Coordinate*  of  Rotating  Helical 
Segment 

radius),  and  ;  (surface  normal).  For  com¬ 
parisons  with  measured  data,  the  cylindrical 
coordinates  (r,  h,  and  z)  are  used  where  r  is 
the  dimensional  radius  varying  in  the  radial 
direction,  u  =  300*;/180  radians  varying  in  the 
r.hordwise  direction,  and  z  varies  along  the 
surface  normal  as  ;  ■ 

Lak shmi narayana ,  et  al.  (1972)  tested  the 
blade  in  a  housing  with  a  9.2(1  cm  ( 0 . 08 in.) 
clearance  between  the  blade  tip  and  the  wall. 

The  free- stream  onset  flow  was  zero  and  the 
rotational  speed  of  the  blade,  denoted  :i , 
remained  constant  at  47  rad/sec  (450  rev/min). 
The  kinematic  viscosity  j  equaled  1.49  x  10"5 
m^/s  (160  x  10"6ft^/sec)  which  corresponds  to 
air  at  20nC  (68°F).  These  conditions,  which 
yield  a  Reynolds  number  based  on  tip  radius  and 
tip  rotational  speed  of  7  x  10’,  are  used  in 
the  evaluation  of  the  present  method.  In  addi¬ 
tion,  present  calculations  begin  at  the  blade 
surface  location  h  =  0.016  radians  and 
xR  =  r/R  =  0.492. 

Initially,  a  boundary-layer  calculation 
of  the  rotating  blade  was  made  with  transition 
set  to  1.0  radian  to  allow  a  comparison  with 
the  laminar  flow  predictions  of  Banks  and  Gadd 
(1962)  on  a  flat  sector  and  of  Morris  (1981)  on 
a  helical  blade.  The  comparisons  with  the  pre¬ 
dictions  of  Morris  (1981)  are  valid  at  large 
radii;  the  twisted  blade  used  by  Morris  has  not 
yet  been  evaluated  with  the  present  method. 
Nonetheless,  the  large  radii  comparison  is 
important  for  validating  the  radial  skin  fric¬ 
tion  coefficient.  Figures  4  a  through  c  show 
comparisons  with  Morris's  (1981)  calculations  at 
r/R  x  0.95.  Figures  4a  and  4b  compare  the 
streamwise  and  radial  skin  friction  coef¬ 
ficients 
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In  these  definitions,  sc  and  \r  are  the  shear 
stresses  in  the  chordwise  and  radial  direc¬ 
tions,  respectively,  »  is  the  fluid  density, 
and  uref  is  the  reference  velocity.  The  solid 
lines  in  rigures  la  and  4b  denote  the  Morris 
computation  and  the  circles  denote  the  present 
calculations.  The  agreement  between  the 
two  prediction  methods  is  excellent. 
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The  solution  of  the  first-order  equations 
is  obtained  using  centered  finite  differences 
0 n  a  net  rube  and  Newton's  method.  Details  of 
• ” ' s  procedure  may  be  found  in  Cebeci  and 
Bradshaw  (19/7). 

The  calculation  of  the  blade  boundary- 
layer  proceeds  as  follows.  Solution  marching 
begins  at  the  hub/ 1 eadi ng-edge  intersection  and 
mores  outward  along  increasing  radii  to  the 
tip.  Calculations  continue  in  this  manner  for 
successively  increasing  chordwise  stations 
until,  finally,  the  blade  trailing  edge  is 
reached.  Calculations  are  made  over  the  entire 
blade  surface. 

figure  2  illustrates  the  calculation  pro¬ 
cedure  in  more  detail.  Calculations  at  the 


leading  edge  begin  as  laminar  at  the  node 
denoted  (A)  In  Figure  2  and  march  radially  out¬ 
ward.  After  completing  the  leading  edge,  the 
solution  is  obtained  at  the  location  of  the 
nest  chord  and  first  outward  radius  from  the 
hub,  node  (B)  in  Figure  2,  using  a  modified 
characteristic  box  method  (Bradshaw,  et  al. 
(1981)).  The  characteristic  box,  which  estima¬ 
tes  a  solution  based  on  the  results  at  the  pre¬ 
vious  chordwise  position  lying  along  the  sane 
streamline,  is  traditional ly  applied  only  in 
regtons  of  reverse  crossflow.  However,  this 
method  has  been  found  to  be  extremely  benefi¬ 
cial  for  propeller  blade  application  near  the 
hub  in  that  it  eliminates  a  complex  solution  of 
hub  initial  conditions. 

Once  a  solution  is  obtained  at  location 
(B),  the  solution  is  set  at  the  current  chord- 
wise  hub  position  using  the  initial  conditions 
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Computations  continue  at  increasing  radial 
positions  up  to  the  tip  using  the  general  box 
in  regions  without  reverse  crossflow  and  the 
characteristic  hox  in  regions  with  reverse 
crossflow.  The  tip  node  is  always  computed 
using  the  general  box.  This  general  solution 
pattern  is  continued  along  the  entire  blade. 

A  transition  computation  is  not  incor¬ 
porated  in  the  present  houndary-i ayer  scheme 
and  the  location  of  transition  to  turbulent 
flow  is  an  input  parameter  supplied  by  the 
user.  The  location  of  transition  nay  be  esti¬ 
mated  by  experimental  results,  such  as  flow 
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Fig  2  Computational  Procedure 
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For  the  special  case  of  «e  «  0,  as  In  the 
rotating  disk  solution  reduces  the  simi¬ 
larity  relations  reduce  to 
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This  hub  similarity  condition  is  exact  for  the 
cases  of  laminar  flow  on  a  rotating  segment  and 
a  two-dimensional  swept  wing. 

By  substituting  the  general  similarity 
formulations  In  Equation  19  Into  the  momentum 
Equations  13  and  14,  the  derivative  quantities 
3w/3n  are  eliminated  and  marching  can  proceed 
in  the  E-direction.  Details  of  the  numerical 
procedure  used  in  the  computer  code  at  the  hub 
are  given  in  Section  6  in  the  numerical  analy¬ 
sis.  The  procedure  is  self-starting  and 
stable. 

Ia_M.BUI L E_NCE_  MODEL 

The  governing  flow  equations  presented  in 
Section  4  contain  more  unknown  quantities  than 
equations.  A  standard  procedure  in  boundary- 
layer  solutions  is  the  introduction  of  a  tur¬ 
bulence  model  to  approximate  the  Reynolds 
stress  terms  -u  v  and  -v'w'.  Various  tech¬ 
niques  have  been  developed  to  handle  turbulence 
modeling,  all  involving  some  degree  of  correla¬ 
tion  with  experimental  data.  The  zero-equation 
model  is  the  simplest  approach  and  does  not 
require  the  solution  of  any  additional  dif¬ 
ferential  equations.  This  method  assumes  the 
algebraic  specification  of  both  length  scale 
and  turbulence  energy  as  explained  and  used 
by  Tebeci  and  Smith  (1974).  The  one-equation 
model  solves  an  equation  for  the  turbulence 
energy  but  uses  an  algebraic  specification  for 
the  length  scale.  Although  this  formulation 
has  been  used  by  Bradshaw,  et  al .  (1976)  for 
thin  shear-layers  with  considerable  success, 
its  usage  is  not  wide-spread.  Two  equation 
turbulence  models,  particularly  the  K-e  model, 
are  experiencing  increased  popularity.  These 
models  are  general,  but  require  the  solution  of 
two  differential  equations. 

For  the  present  work,  the  simple  zero- 
equation  model  of  Cebeci  and  Smith  (1974)  has 
been  adopted.  An  eddy-viscosity  factor  . m  is 
defined  to  relate  the  Reynolds  stresses  to  the 
mean  velocity  profiles  by 


-u  v'  -  cm  yi  and  -v'w'  =  cm  (21) 


The  theoretical  boundary  layer  is  divided 
into  an  Inner  region  and  an  outer  region  with  a 
separate  equation  defining  em  In  each  region. 
The  Inner- region  formulation  Is  applied  from 
the  wall  surface  to  the  location  In  the  boun¬ 
dary  layer  where  both  inner  and  outer  equations 
predict  the  same  eddy  viscosity.  From  this 
location  to  the  edge  of  the  boundary  layer,  the 
outer  wall  formulation  Is  applied. 

Eddy  viscosity  In  the  inner  region  Is 
defined  as 
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where  L_is  a  mixing  Jength  approximation  equal 
to  0.4c  l-exp(l-c/A)  ,  A  is  a  damping  factor 


equal  to  26 ,  and  Ttw  is  the  turbulent 

1.  e  ) 


wall  shear  stress  equal  to  u  +7i5\ 

■’-.,'w  ?  c  4 

+2cosa|.au'  ^jwn  )£  .  in  the  outer  region,  the 
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eddy  viscosity  becomes 

(Em)  =  0.0168  I  ■  (u^p  -  U{ ) dc  1  i^r 
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where  u^e  equals  (ue?  *  we?  *  2uewecosa))£, 
u*  equals  (u?  ♦  ♦  2uwcosu)‘£  and  is  an 

intermittency  factor  to  account  for  the  tran¬ 
sition  region  between  viscous  and  turbulent 
flow. 


6.  NUMER_ICAL_  ANALYSIS 

The  governing  three-dimensional  boundary- 
layer  equations  for  propeller  blades  are 
solved  numerically  using  the  Keller  (1910)  bo* 
method.  This  solution  technique  may  be  divided 
into  four  steps.  Initially,  the  governing 
equations  are  written  as  a  system  of  first 
order  equations  by  the  introduction  of  trans¬ 
formed  variables.  The  first  order  equations 
are  then  written  in  finite  difference  form 
using  central  differences.  Newton's  method  is 
applied  to  linearize  the  difference  equations 
and,  finally,  the  linear  system  is  solved  by 
the  hloclt  triangular  elimination  method. 

New  variables  u,  v,  w,  t,  and  ',  defined  as 

i '-  v  =  f  " 


w '  *  t.  =  g  " 

‘3  =  mju  *■  mq,w  *  P]i)  i-  1117  24  ) 


are  introduced  to  reduce  the  •*.-  and  ■  -momentum 
equations  to  first  order.  The  new  system  is 

u  =  v 

w  =  t 

f  bv  I  -  mou^  -  intjiiw  -  m^w<?  *  in j  j 
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The  new  terras  1*13  through  1*17  account  for  the 
blade  rotation  and  are  defined  as 

i*13  =  2ilsi  njCOSaCSCa 
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The  boundary  conditions  for  the  transformed 
equations  are 

f  *  f'=  g  ■  g'  *  0  at  0  =  0 
f'=  1,  9'=  Webref  a*  0  =  °- 


4.  INITIAL  CONDITIONS 

The  propeller  boundary-layer  solution  is 
obtained  by  marching  In  the  radial  direction 
for  chordwlse  locations  Increasing  from  the 
leading  edge  to  the  trailing  edge.  In  this 
manner,  the  computation  covers  the  entire  blade 
surface.  The  solution  procedure  requires  ini¬ 
tial  velocity  profiles  along  two  intersecting 
planes.  These  planes  are  chosen  to  be  the  pro¬ 
peller  leading  edge  and  hub.  Although  the 
marching  begins  at  the  radial  location  termed 
the  hub,  the  actual  flow  at  hub/blade  intersec¬ 
tion  Is  too  complex  to  be  represented  by  the 
boundary-layer  analysis.  For  this  paper,  the 
term  hub  defines  a  small  radial  distance  out¬ 
board  of  the  actual  hub/blade  intersection 
where  the  boundary-layer  equations  apply. 

Since  the  Initial  conditions  are  not,  in 
qeneral ,  known  quantities  for  a  given  blade 
geometry,  assumptions  become  necessary  to  begin 
the  computation.  The  remainder  of  this  section 
describes  the  initial  conditions  adopted  in  the 
present  propeller  boundary-layer  calculation. 


For  the  rotating  helical  blade.  Banks  and 
Gadd  (1962)  theoretically  show  that  the 
leading-edge  similarity  function  for  the 
streamwlse  velocity  u  satisfies  the  Blaslus 
equation.  This  analysis  can  be  extended  for 
application  to  the  propeller  blade  leading-edge 
using  the  equations  for  the  similarity  solution 
of  wedge  flow.  The  leading-edge  solution  of 
this  paper  Is  similar  to  that  adopted  by  Itoh, 
et  al .  (1984)  In  which  the  Karman-Pohl hausen 

profile  is  used.  However,  the  current  method 
computes  the  crossflow  velocity  component  w  at 
the  leading  edge  from  the  n-momentum  equation. 
Itoh,  et  al .  (1984)  set  these  profiles  to  zero. 
With  the  present  analysis,  the  governing 
equations  at  the  leading  edge  are 

r, -momentum:  (bf")'+  miff"-  i*2(f')^  +  mu  =  0 

n-momentum:  (bg")'+  r*ifg"+  *132  =  0 

with  the  same  coefficient  terms  b,  mj,  mj,  mu 
and  1*12  defined  in  equations  15  through  17. 
These  equations  compute  two-dimensional  stagna¬ 
tion  flow,  wedge  flow,  or  flat  plate  flow  at 
zero  Incidence  depending  upon  the  specified 
external  potential  flow  velocities  at  the 
leading  edge.  This  solution  has  proven  to  be 
both  stable  and  smooth  In  the  boundary-layer 
calculation  at  the  leading  edge  and  does  not 
require  averaging  or  Iterating. 

For  nonrotating  three-dimensional  calcula¬ 
tions,  the  second  Initial  condition  is  spe¬ 
cified  at  locations  where  the  crossflow 
velocity  w  equals  zero  everywhere  Inside  the 
boundary  layer.  For  the  ship  hull  (Cebeci ,  et 
al  .  119781)  and  the  arbitrary  wing  (Cebeci  et 
al .  (1977  )),  the  locations  of  zero  crossflow 
are  the  ship  keel  and  the  wing  root,  respec¬ 
tively.  For  these  regions,  the  n-momentum 
equation  is  Identically  equal  to  zero  and  the 

d  W 

unknown  variables  become  u,  v,  and  —  -  .  That 

an 

is,  for  w-0.  Equation  14  becomes  0=0. 

To  obtain  a  solution,  a  new  equation  is  derived 
by  differentiating  the  singular  n-momentum 
equation  with  respect  to  n.  This  equation  is 
called  the  attachment  line  equation. 

For  a  propeller  blade,  the  solution  w:0 
everywhere  inside  the  boundary  layer  is  valid 
only  for  the  blade  radius  equal  to  zero  and, 
therefore,  should  not  be  used  to  begin  the  com¬ 
putation  for  arbitrary  nonzero  radii.  Itoh,  et 
al  .  (1984)  apply  the  attachment  line  solution 
to  the  huh  of  the  axial  flow  pump  blade  hut 
find  it  necessary  to  use  an  Iterative  method  to 
obtain  a  reasonable  solution  at  the  hub.  In 
this  paper,  the  method  of  a  similarity  solution 
is  adopted  at  the  hub  as  well  as  at  the  leading 
edge.  The  similarity  relations  used  at  the 
present  time  are: 
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